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Abstract 

Based on a selfconsistent theory of localization we study the electron transport properties of a 
disordered system in the framework of the Anderson model on a Bethe lattice. In the calculation of 
the dc conductivity we separately discuss the two contributions to the current-current correlation 
function dominating its behaviour for small and large disorder. The resulting conductivity abruptly 
vanishes at a critical disorder strength marking the localization transition. 
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Disorder strongly affects the motion of an electron, and even may fully suppress electron 
transport, as first revealed by P. W. Anderson^. The natural quantity to characterize the 
effect of disorder is the electric conductivity, whose calculation for a disordered interacting 
system should be understood as a primary goal in the field of localization physics. As a first 
preparative step, we consider the dc conductivity for a single (noninteracting) electron in 
a crystal with compositional disorder. Our work is directed towards an approach which is 
applicable for all values of disorder. 

The generic model in the field of localization is the Anderson model 

ff=X>i*><ii-j]Cii>oi • c 1 ) 

Here — J ^2 * s a tight-binding term with hopping matrix element J on a certain lat- 
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tice and X/ e *b")(*1 describes random local potentials. The e« are assumed to be identically 
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independently distributed (i.i.d.) variables with distribution p(ej) = (1/7)0(7/2 — |ej|). 7 
specifies the strength of the disorder. 

For noninteracting particles at T = 0, the Fermi function f(u) in the Kubo formula 
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for the dc conductivity becomes a 5-function, and o"d c is solely determined by the disorder 
averaged current- current correlation function 

XjjM = ^2 ' £')( ImG i,j+<5'M lmG j4+s (uj)} av (3) 
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at the Fermi energy Ep (for interacting particles, this corresponds to a rigid band approx- 
imation). Gij(u) denotes the retarted single-particle Green function. Below, we assume 
Ep = (half filled band). 

To avoid the need to sum over all lattice sites in eq. (J3J), Girvin and Jonson^ suggested 
to split this expression in two parts, 

X jj(u) = P!(u)A(u), (4) 

with the disorder averaged pair correlation function 

Pi = (1m Gu(uj) ImGj^w) - ImG^u) ImG ii (u;)) av (5) 
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for adjacent lattice sites i, j, and a correction A(u) accounting for the long range correlations. 
Note that in the regime of localized states A(u;) has no definite value since both Pi (a?) and 
Xjj vanish. Close to the localization transition xjj is essentially determined by the pair 
correlation function, and A(uj) is nearly constant. For small disorder, on the other hand, 
A (a;) diverges. To cover the full range of disorder, we must therefore compute both P\{oS) 
and A(to>) to a good approximation. 

To calculate the pair correlation function Pi(cu), which captures the behaviour near the 
localization transition, we employ the self-consistent theory of localization developed by 
Abou-Chacra, Anderson, Thouless (AAT) 3 . This theory is based on a renormalized pertur- 
bation expansion 4 which sets up a closed set of recursion relations for local Green functions 
on the Bethe lattice. These recursion relations can be interpreted as a stochastic selfcon- 
sistency equation for the Green function Gjhu), which corresponds to the lattice with site 
i removed (and appears in the second step of the renormalized perturbation expansion). 
Solving this stochastic equation by a Monte-Carlo procedure a sample for G^{u) is con- 
structed from which the respective distribution can be calculated. On the Bethe lattice the 
full Green function as well as all nondiagonal Green functions may be expressed in terms of 
local Green functions G^{uj) as 

G«M = —T-, n , Ga(u) = JGdu)GfAu) , (6) 

with j nearest neighbour i. Furthermore Pi(&>) can be calculated from the distribution of 
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Note that, since the Bethe lattice has no closed loops, G^{u) and G^{oS) are i.i.d. random 
variables. 

In figure ^ we depict, for the Bethe lattice with connectivity K = 2 and for bandwidth 
W — 1, the pair correlation function Pi (a;) at the band center (u = 0). For not too 
large disorder, P\{oj) is finite, but abruptly vanishes for 7 > 2.9, indicating the transition 
from extended to localized states. As figure moreover shows the behaviour of the pair 
correlation strongly depends on the imaginary part 77 in the energy argument u + ir] of the 
Green functions. To detect the localization transition it is necessary to perform the limit 
77 — > numerically (for details on this, see 5 ). 
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FIG. 1: Pair correlation Pi(to) at the band 
center (u> = 0). 7 is measured in units of the 
bandwidth. 77 denotes the imaginary part of 
the energy argument in the Green functions. 



FIG. 2: Comparison of typical DOS and pair 
correlations for two values of 7. The solid 
line shows the pair correlations Pi(cj), the 
dashed line the typDOS. Results are rescaled 
to their value at u = 0. 



Figure El shows a comparison between the pair correlation and the typical density of states 
(typDOS) 



N typ (u) =exp((lniVj; ) H) 
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with Njj(uj) 
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which has been suggested as an order parameter for localization^. Remarkably enough, the 
pair correlation follows the behaviour of the typDOS over the full range of energy, implying 
a close relation between these two quantities. In this sense, the typDOS might be itself 
understood kind of transport quantity. 

In contrast to Pi(uj), the correction A(u) contains contributions from all terms in eq. Q. 
Since directions on the Bethe lattice are ill-defined a direct calculation of these terms within 
AAT is not possible. A natural suggestion to overcome this obstacle is to use the coherent 
potential approximation (CPA). While for small disorder the CPA surely produces a senseful 
result the situation for large disorder is less clear. To check whether a CPA calculation of 
A(o>) close to the localization transition is reasonable we have to estimate the qualitative 
behaviour of A (a;) within AAT. Running along a single non-retracing path in the Bethe 
lattice, all terms in eq. Q along this path have the form 



P n (u) = (ImG iii+ri _i(w) lmG i+1)i+n (uj) - ha.G iti+n (u) ImG , i+ i ii+n _i(^)) £ 
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similarly to the one-dimensional chain (of course, for n = 1 we get back the pair correlation 



function). On the Bethe lattice (K + l)K n 1 paths of length n originate from a single site, 
so that the magnitude of A(u;) can be estimated as^ 

AM-^+ir- 1 ^ . (io) 

Without disorder (7 = 0) each term decays as P n {uj) ~ K~ n , and the sum diverges. With 
disorder, the P n (co) decay faster, resulting in a finite value for A(u). Figure El demonstrates 
that A(u) is non-zero and varies only slowly in the vicinity of the localization transition 
(7 C w 2.9). Besides confirming the original expectation, this indicates the possibility to 
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FIG. 3: A(w)-correction to the pair correlation functions at the band center (u = 0). The solid 
line is the AAT-estimate according to eq. (jlUj) . the dashed line shows the CPA-result. The vertical 
dotted lines indicate the position of the localization transition at 7 C ~ 2.9. 

employ the CPA for the calculation of A(z). From the CPA current-current correlation 
function Xjj A (u;)£ the correction A(cu) can be approximated as A CPA = Xjj A /Pi PA (inset in 
figurelHJ). That A CPA (u;) is non-zero even above the localization transition makes no problems 
in our approach: Since the pair correlation P\{oj) is zero, xjj^) is as weu ) indicating 
localization of states. 
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Now we are in the position to put pieces together: If we calculate the pair correlation 
function P\{oS) within AAT, and the correction A(u;) within CPA, the current-current corre- 
lation function xjj(cj) is expected to be obtained in good approximation. Although it is not 
necessary to split the current-current correlation function in the proposed way, it seems quite 
reasonable. First, we have identified the two contributions to the current- current correla- 
tion function which dominate its behaviour in the regimes of small respective large disorder. 
Second, since both contributions show different behaviour we can apply different methods 
for their calculation. 

The pair correlation function P\{oS) rules the behaviour of xjjI^) f° r large disorder. 
While it is finite for all values of disorder, it is critical at the localization transition. Indeed 
the full effect of localization is contained in P\{uS) but not in A(z). Mean field treatments 
like the CPA, merely focusing on average values, are not sufficient for P\{u). The ap- 
parent relationship between the typical DOS and the pair correlation function makes this 
point especially clear: The pair correlation function does not agree with the averaged DOS 
("arithmetic mean"), but with the typDOS ("geometric mean"). This failure of mean field 
treatments indicates the subtlety of the problem at hand. In constrast the selfconsistent 
theory of localization by Abou-Chacra et. al.— can be conveniently applied, and allows for 
a not too complicated calculation of the pair correlations. The correction A(u) is in some 
sense opposite to the pair correlation function. While it is noncritical at the localization 
transition, it diverges for zero disorder. If one is only interested in the behaviour close to the 
localization transition one might completely forget about A (a;). But to correctly describe the 
regime of small disorder A (a;) has to be taken into account. Nevertheless, since A (a;) shows 
no critical behaviour, there is no need to go beyond a mean field treatment, e.g. provided 
by the CPA. As additional benefit, this avoids all possible problems arising from the lack of 
well defined directions on the Bethe lattice. 

In conclusion, the proposed splitting of the current-current correlation function seems to 
be a successful first approximation. Since the pair correlation P\{oj) is a transport quantity 
it is a natural localization criterion, whose calculation allows for a precise determination of 
mobility edges and the critical disorder for the localization transition. The relation between 
P\{ui) and the typDOS shows that the latter one can indeed be used as a localization 
criterion as in^^. The A(u;)-correction, on the other hand, is of no importance to detect 
the localization transition, but accounts for the divergence of the current-current correlation 
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function as 7 — > 0. The final result for Xjj(u;) correctly interpolates between the two limiting 
cases 7 — > and 7 — > y c , where it matches the exact result. 

Note that the AAT and the CPA can be combined to some extent with a treatment 
of interaction processes (statDMFA, DMFA/DCPA). Adopting the ideas presented here to 
these extensions might hint at some possible direction for future attempts to deal with 
electron transport in a disordered interacting system. Especially the selfconsistent theory of 
localization^, which can be straightforwardly combined with a local treatment of interaction 
(statDMFA 5 ' 8,9 ), is a promising candidate for these studies. Work on this subject is in 
progress. 

One of the authors (A. A) acknowledges support from the European Graduate School, 
Bayreuth. 
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